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Abstract. There are no known exact formulas for the valuation of a num- 
ber of exotic options, and this is particularly true for options under discrete 
monitoring and for American style options. Therefore, one usually recourses 
to a Monte Carlo Simulation approach, amongst other numerical methods, to 
estimate the value of these options. The problem which then arises with this 
method is one of variance reduction. Control variates are often used, and we 
present some results concerning these control variables, for the valuation of 
Asian and lookback options. An inequality on functions of correlations useful 
for comparing estimators in variance reduction procedures is also provided. 



1. The problem of exotic options valuations 

There are two main types of options. Standard European or American put and 
call options are referred to as vanilla options. These are exchange traded options 
and their payoffs depend on the price of the underlying asset at exercise time 
only. Contrary to these, there are derivatives with more complicated payoffs and 
which are referred to as exotic options. Most exotic options trade in over-the- 
counter markets and have been designed to meet the particular needs of customers 
(corporations). 

One important class of exotic options consists of path-dependent (or history- 
dependent) derivatives. They have a payoff that depends not just on the final 
price of the underlying asset, but also on the path followed by this price. Asian 
options and lookback options are typical examples of path-dependent options. In 
some instances, numerical methods are the only means available to value these 
exotic options, and in addition to Monte Carlo methods, other commonly used 
numerical methods include the finite difference and the finite element methods, the 
trinomial tree method, and some times the Markov Chain method (see [7, 8, 9, 
10, 11] and the references therein). Until the last decade, Monte Carlo method 
has been considered as costly and unreliable, but based on innovative techniques 
which are now a topic of current research, they've yielded more promising results. 
In this context, Barraquand and Martineau [10], amongst others, obtained some 
interesting results. 

Consider an option on a given security having n days to maturity. Let Sd(j) 
denote the closing price of the security at the end of day j ( for j — 1, . . . ,n). 
Assume also that the risk free interest rate is r per year, and denote by N the 
number of trading days in a year. With these notation, the payoff from an asian 
option is either given by 
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where the strike price X is taken as the average of all prices up to the expiry date, 
or 



when the underlying's terminal price is taken as the average of all prices up to 
maturity. The payoff from a lookback option is given by 



Barrier options are another type of frequently traded path-dependent options. 
These are options for which a barrier value v, which may be a one- or two-dimensional 
vector, is specified, and the option becomes alive or cease to exist according to the 
path followed by the price w.r.t. the barrier value. When a barrier option is alive 
at maturity, its payoff is that of the equivalent vanilla option, i.e. that with the 
same strike price K, and number of days n to maturity. 

The main problem with using Monte Carlo Simulation to value path-dependent 
derivatives is that the computation time necessary to achieve a reasonable level of 
accuracy can be excessively high. Control variate techniques are some of the tools 
available for implementing the necessary variance reduction procedures, which can 
lead to dramatic savings in computation time. 

Our focus in this paper will be on the valuation of Asian and lookback options. 
In section 2 we discuss Monte Carlo valuation simulations for these options and 
explain how controlled variates techniques can be gainfully employed to improve 
the accuracy of the valuation. In section 3, we present some results concerning the 
control variables for an improvement of the variance reduction procedure, and an 
inequality on functions of correlation coefficients useful for comparing estimators 
in variance reduction techniques is also provided. 

2. Monte carlo simulations and variance reduction procedures 

Suppose that the price S(t) of a given security follows a risk neutral geometric 
Brownian motion with constant parameters /i and a 2 . This means in particular that 
for all y > 0, and t > 0, ln(S(t + y)/S(y)) has a normal distribution with mean /it 
and standard deviation a\ft, and that fi = r — a 2 /2, where r represent the risk-free 
interest rate. 

We denote as usual by 5(0) the initial price of the security, and by Sd(j) the 
price of the security at the end of day j ( for j = 1, . . . , n), where n is the number 
of days to maturity. Set 




(1.2) 




(1.3) 
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The random variables X(l), . . . , X(n) are independent and identically distributed 

J± si' 

N ' N ■ 



with mean fi/N and variance a 2 /N. That is, X(i) ~ 3>(-j£, ^r). A straightforward 



calculation using (2.1) shows that 

S d (i) = S d (i - l)e x ^ (2.2) 
= S(0)e x W + - +x ® (2.3) 

Values of -X'(l), • • • , X(n) can be generated by a computer and Eq. (2.2)) or (2.3) 
can be used to sample a random path for the price S of the security. 
The valuation simulation for an asian option whose payoff is 

Y = e ~ rn / N ^2 - K 

as in Eq. (1.2) can be implemented as follows 

(1) Generate n values for the normal random variables X(l), . . . , X(n). 

(2) Calculate the end-of-day prices S d (i) = S(0)e x ^+- +x( - n \ for i = 1, . . . , n. 

(3) Calculate the payoff of the derivative Yj = e - rn / N (j2 - K^j , for 
the j-th simulation run. 

(4) Compute the average mean Y — ji^2i'Yj, where R is the total number 
of simulation runs, i.e. the number of times that steps (1) through (3) arc 
performed. 

The value of Y represents an estimate of the exact cost Y of the derivative. We 
have 



E[y] = E[y] (2.4) 

and 

Var(?) = Vsx(Y)/R. (2.5) 

The last equality, Eq. (2.5), shows that the accuracy of Y improves with the 
number of simulations. 

The simulation procedure is similar for lookback and for most of the path- 
dependent derivatives. Monte carlo simulation presents some advantages. They 
can be used to value path-dependent derivatives as well as those derivatives whose 
payoff depend only on the final value S of the security. They can also be used to 
value barrier options type securities, whose payoffs may occur several times before 
maturity. More discussions on the advantages of Monte Carlo simulation techniques 
can be found in [1, P. 408]. 

The drawbacks with Monte Carlo simulation is that the amount of time nec- 
essary to achieve a reasonable accuracy can be unacceptably high. As remedial 
measures, a number of variance reduction techniques are available, and they can 
lead to dramatic savings in computation time. These include the antithetic vari- 
able technique, the importance sampling, the stratified sampling, and the control 
variate technique. All these procedures and many others are described in [1, 2, 3]. 
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We now move on to explain how control variate can be used to improve the 
simulation valuations of derivatives. Suppose that in a general setup we want to 
estimate 



M = E[Y] (2.6) 

by simulation. Suppose furthermore that in the process of generating the random 
variables that determine Y we also find another random variable V whose mean 
value is known to be /j, v = E[V]. Then rather than using the value of Y as an 
estimator, we can replace it with another value 



W = Y + c{V-fi v ), (2.7) 

where c is a constant to be specified. The value of W is an estimate for M since 
E[W] =E[Y] = M. We have 

Var(W) = Var(y ) + c 2 Var(V) + 2c Cov(Y", V) (2.8) 

The best estimator is obtained by choosing c to be the value that makes Y&v{W) 
as small as possible. Elementary calculus shows that the value of c that minimizes 
Var(W) is given by 

,*- 9^X1 (29) 
~ Var(F) 

substituting this value back into Eq. (2.8) yields 

v^-W,-^. 

That is, 

(2,0, 

= l-r 2 , (2.11) 

where the quantity r = Corr(V, V) is the correlation coefficient between Y and 
V. The number r 2 is often called the coefficient of determination between Y and V. 
Eq. (2.10) shows that the variance reduction obtained with the control variate V 
is completely determined by Corr(Y, V) . More precisely, the reduction obtained is 
100% Corr 2 (F, V) and consequently, the higher the correlation between Y and V, 
the greater the variance reduction obtained will be. 

2.1. Control variates. Suppose that the payoff Y from an option on a given 
security is of any of the forms given by the Eqs. (1.1)-(1.3). With the notation of 
the previous section, set 



V = J2^X(i) (2.12) 

i=l 

where at are some weights to be determined, and take 
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W = Y + c(V - E(V)) 



as the new estimator for Y. 
We have 



so that 



where 



/ 2 \ 

E[F]=^a i ETO]=^|(r-^), 



i=l 



VK = y + c 
= r + £/* 



5>*(i)-£(r-^-) 



X(»)-(r-^/JV 



(2.13) 



Pi = cai. (2.14) 

On the other hand, we have 

n n 

Var(VK) = Var(F) + 01 Var(X(i)) + 2 ^ ft Cov(y, 
and it readily follows that the values of ft that minimize W are given by 

Cov(y,x(i)) 

* " Var(X(i)) (2 ' 15) 

In terms of these optimal values of the ft, the smallest possible value for Var(iy) 
is given by 



That is 



(2.16) 



= l-^Corr 2 (r,XW) 



(2.17) 



Consequently, 

Var(W) 
Var(F) 

We have thus proven the following result 

Theorem 1. Let V = ^ i aiX(i), for some arbitrary weights on and fori = 
1, . . . , n. Let W = Y + c(V — E[V]) be an estimator of Y for some constant c. 
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(a) : The optimal variance reduction is achieved with this control variable V 
for the values of on and c such that 

Cov(Y,X(i)) 
Ca * = - Var(X(z)) ' (• = 1 '-.») 

(b) : The optimal variance reduction thus obtained is given by 

Remark 2. 

(1) It follows from part (a) of the theorem that a control variate of the form 
V — J2i a iX(i) gi ves r i se to an optimal reduction if and only if on = 
~^ C °Var(x(i)' > f° r alH = 1, . . . , n and for a nonzero constant A. 

(2) This theorem implies that for all a\, . . . , a n and for every 
constant c 



where 



W = Y + c 



J2<XiX(i)- ai (r-o- 2 /2) 



In [4, P. 138], the problem of finding the values of «j ( for i = 1, . . . , n) for the 
best control variate of the form V = J2i a iX(i), and alternatively, the problem of 
finding the values of c\ , . . . , c n for the best estimator of Y of the form 

r 2 , 



i=l 

are posed. Theorem 1 gives an answer to this question, by determining the values 
of on and a and by indicating precisely the optimal variance reduction that can be 
achieved. 

Theorem 3. Let Y be a random variable and let X\, . . . ,X n be n independent 
random variables. Then for all numbers ol\,. . . ,a n 



Corr 2 (Y,J2^x)j < ^ Corr 2 (y, X t ). 



Proof. Let V = ^ and let W = Y + c(V - E[V]) be an estimator of Y for 
some constant c. For any fixed values of the Qj the optimal variance reduction given 
by Eq. (2.10) is 

^"'-^ 

This is not smaller than the optimal variance reduction obtained for all possible 
values of the on and c given by Eq. (2.17). Consequently, 1 — J2 i Corr 2 (F, Xi) < 
1 - Corr 2 (y,E 4a4 ^); that is, Corr 2 (F, ^ ^Xi) < E 2 Corr 2 (F, X t ). □ 

The result stipulated by this theorem is certainly very important, but it is likely 
to be unknown. Indeed in statistics books that present the most comprehensive 
material on the topic of correlation coefficient between random variables, there is 
rarely any discussion of a relationship between functions of correlations (see [5, 6]). 
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2.2. Applications. 

(1) If V = J2i a iXi where the a.i are constants independent of Xi, then no 
estimator of Y of the form W = Y + ctiXi — E[V] using V as control 
variable can lead to an optimal variance reduction, since by Theorem 1, the 
on must depend on Xi. 

(2) Let V = Ln(s d (n)/S(0)). By the definition of the X(i), it readily follows 
that V = ^2iX(i). Thus by the preceding remark, taking V as control 
variable cannot give rise to the best estimator. 

3. Conclusion 

We have made use of the variance reduction formula (2.10) obtained for an 
estimator of the payoff Y, with a given control variable V, to derive the form of 
the best possible control variable of the form V = J2iCtiX(i), where the random 
variables X(i) are as in Eq. (2.1). We've also obtained an upper bound for the 
correlation coefficient between an arbitrary random variable Y and an arbitrary 
linear combination of independent random variables. This inequality can be useful 
in comparing estimators in variance reduction procedures. 

A generalization of this type of inequalities between functions of correlation 
coefficients is desirable, to compare for instance an estimator obtained with a control 
variable of the form V = J2iPi^d(i), with that obtained using a control variable 
of the form Y]. ajX(i). In this instance, one would need a general relationship 
between functions of correlation coefficients of the type Corr (Y, 0iSd{i)) and 
Corr 2 (F, £\ Pi H S d(i)/S d (i - 1))). As stated in this paper, no result on such type 
of relationship seems to be available in the literature. Finally, the results obtained 
in this paper clearly apply to any random variable for which variance reduction is 
required, and not only to the specific case of exotic options considered in this paper. 
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